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Abstract 

I present a new improved estimator for the correlation function of 
2D nonlinear sigma models. Numerical tests for the 2D XY model and 
the 2D 0(3) -invariant vector model were performed. For small physical 
volume, i.e. a lattice size small compared to the bulk correlation length, 
a reduction of the statistical error of the finite system correlation length 
by a factor of up to 30 compared to the cluster-improved estimator was 
observed. This improvement allows for a very accurate determination 
of the running coupling proposed by M. Liischer et al. for 2D 0{N)- 
invariant vector models. 
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1 Introduction 



Almost twenty years ago Nightingale pointed out that the correlation length of 
a finite system can be used to trace the renormalization group fiow of a statistical 
model. He demonstrated the validity of his approach at the example of the 2D Ising 
model using its exact solution 0. 

Liischer computed the mass gap (the inverse correlation length) for the 2D 0(3) 
vector model on finite lattices using perturbation theory 0. In ref. 0] Liischer, 
Weisz and Wolff used the mass gap of finite systems to define a running coupling 
for 2D asymptotically free models. The aim of their study was to match the pertur- 
bative, i.e. small distance, regime with the nonperturbative large distance regime, 
using MC simulations for large distances and in the cross-over region. The high 
numerical accuracy that they reached for the step scaling function of this particular 
running coupling of the 2D 0(3)-invariant vector model was based on the use of the 
efficient cluster algorithm ||^ combined with a variance reducing cluster-improved 
estimator [^]. However one has to note that up to now this powerful simulation 
technique only applies to a rather restricted class of models. I will discuss an al- 
ternative improved estimator that does not rely on the cluster algorithm and hence 
can be applied to a larger class of symmetry groups. 

The concept of improved estimators is rather general and by no means restricted 
to the cluster algorithm. The basic idea of improved estimators is to integrate ana- 
lytically or by standard numerical methods those degrees of freedom that are most 
important for the fluctuations of the observable in question. For local fluctuations 
this can be achieved easily. Parisi et al. proposed to improve the estimator of 
the Wilson-loop correlation-function by integrating over the link-variables on the 
Wilson-loops. 

I will show that such ideas also apply to large range correlation functions of 
certain models. The estimator is based on an embedding of a one-dimensional model. 
The degrees of freedom of this embedded model are collective transformations of all 
spins in a given time-slice. 

This paper is organized as follows. In section 2 I deflne the new improved esti- 
mator for two-dimensional 0{N) invariant vector models. In section 3 I summarize 
theoretical predictions for the running coupling of 2D 0(A^)-invariant vector models. 
In section 4 I present numerical results for the 2D 0(2)-invariant vector (XY) model 
and the 2D 0(3)-invariant vector model. The conclusions are given in section 5. 
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2 The improved estimator 



First let us discuss the strategy applied for finding an improved estimator for an 
observable A of a statistical model with field variables 4>. The expectation value of 
A is given by 

_ /£>0exp(-5(0)M(0) 

JD<l>eM-Sm ■ ^ ' 

We split the set of field variables into two parts (f) — (0, $), where $ is chosen such 
that the integrals over $ can be performed exactly for any fixed 0. Let us define 
the improved estimator by 

A^mA^) - - (^)cond " / exp(-^(0, $)) ' 

where {...)cond denotes the conditional, i.e. dependent expectation value. It is easy 
to check that {Aimp) — and also that the variance of the improved estimator 
Aijnp is reduced compared to that of A 

{A') - {Af > (AL,) - {Arrnp? • (3) 

Of course, the remaining non-trivial question is to identify, for a given model and 
observable, the $ such that the integration over $ can be performed with reasonable 
effort and leads to a considerable reduction of the variance. 

Now let me explain the improved correlation function estimator for the 0{N)- 
invariant vector models. The action is given by 

N 

S^-PT.T. t)^o.{x +l,t)+ Sa{x, t)Sa{x, t + 1)) (4) 

x,t a=l 

where Sa{x, t) are the components of a real A'"-component unit vector. The lattice has 
extension T in the time (t-) direction and extension L in spatial (x-) direction. For 
simplicity we set the lattice-spacing a equal to 1 in this section. In time direction we 
employ open boundary conditions while in x-direction periodic boundary conditions 
are employed, i.e. identify x = L + 1 with x = 1 and a; = with x = L. We choose 
left-multiplications of all spins in a time-slice t with the same 0(A^)-matrix R{t) 
as degrees of freedom of the ID embedded model. The improved estimator is then 
constructed by integrating over the R{t) exactly. The conditional, i.e. Sq, dependent, 
action of the embedded models is given by 

N / N \ / ^ \ 

ScondiR) = E E E RaAt>s{x, ^) E ^0.^^ + 1)^7(^> ^ + 1) (5) 
x,t a=l \S=1 ) Yi=\ I 
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Note that contributions from spatial bonds stay constant under the transformation, 
and can hence be ignored. The conditional partition function is given by 

Zcond = j DRex-p{-Scond{R)) , (6) 

where DR = YitdRit) and dR denotes the Haar measure of the 0{N) group. Next 
let us introduce relative rotations X{t,t + 1) of two neighbouring time-slices by 

X(t,t + 1) ^ R-\t)R{t + l) (7) 

All R{t) can be expressed in terms of R{1) and the X{t, t + 1). 

R{t) = R{1)X{1, 2)...X{t -2,t- l)X{t - 1, t) (8) 

Now it becomes clear why free boundary conditions in time direction are used. 
Periodic boundary conditions in time direction would require that the product of all 
X{t,t + 1) on the lattice is equal to the unit matrix. In terms of the X the action 
of the conditional model is given by 

N N 

ScondiX) = - E E E ^-,-y{t, t + l)Qa,-y{t, t + 1) (9) 
t a=l 7=1 

where Q is given by 

Qa,j{t, t + 1) ^ P^Sa{x, t)s^{x, t + 1) . (10) 

X 

The properties of the Haar measure allow also to rewrite the measure in terms of 
the relative rotations X 

Zcond = / dR{l) J DXeM-Scond{X)) , (11) 

where the integration over R{1) just gives a constant factor that hence can be 
ignored. The conditional partition function factorizes in two-slice partition functions 

t=T-l 

Zcond = n ^(^,^+1) (12) 
t=l 



with 



j(i,i + l)- /dXexp (-EE^">7(^>^+l)Qa,7(^,^ + l)| ; (13) 
•' \ a=l7=l / 
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where dX denotes the Haar measure of the 0{N) group. 

Next let us turn to the measurement of time-shce correlation functions. The 
time-slice correlation function is defined by 

G{T)=J2Sa{t)Sa{t + T), (14) 
a 

where 

Sa{t)^J2^c.{x,t) . (15) 

X 

The improved estimator of the time-slice correlation function is given by 

Gimp(r) = {G{r)),ond = E(E RaAt)Ss{t) E RaAt + r)S^(t + r))^^^ , (16) 

as 7 

where 

^cond 

is the conditional expectation value of an observable A{R). Again we can re-express 
things in terms of relative rotations 



Gimpir) = J212{Sa{t)X^At^t + T)S^{t + T))^^ (18) 

a 7 

= j:T.Sa{t){x^,,{t,t + r))^^^S,{t + T) , 



a 7 



where X{t,t + t) — R ^{t)R{t + r). Using the laws of matrix calculus and the 
factorization of the conditional partition function we obtain 



Hence the evaluation of the conditional expectation values is reduced to the evalu- 
ation of the low dimensional integral 

X(ff^u\ = I dX expjj:^^^ X^^^Q^At,t + 1))X 

2.1 Performing the integrals for N=2 and 3. 

For N — 2 this integral can be solved in closed form. One obtains 

X„=X.. = ^M*i±«H (21) 

Io{k) k 



and 

V _ Y _ -^l(^) - Q12 

A2I — — A12 — , . , [^^ 

where 



1^ = ^{Qu + Q22? + (Q21 - gi2)2 (23) 

and /q and /i are modified Bessel-functions. 

For = 3 I could not find an exact solution for the integrals. For practical 
reasons I restricted the integration to the 5*0(3) subgroup of 0(3). 

Using the property 

j dXfiX) = J dXfiUXV) (24) 

of the Haar measure, where U and V are elements of the group manifold, the inte- 
grand 

J rfXexp(Tr XQ^) = J rfXexp(Tr UXVQ'^) = J rfXexp(Tr XVQ^U) (25) 

can be transformed such that the source Q = VQ^U becomes diagonal. 

Now I parametrized the 5*0(3) using the Euler-angles. One of the 3 integrations 
can be performed analytically. The remaining 2D integral I solved numerically. 

If all three sums qi = Qn + Q22, q2 = Qii + Q33 and qs = Q22 + Q33 are positive 
the integrand takes its maximal value for X equals the unit matrix. I performed 
a saddle-point approximation around this maximum up to order q^^. Comparing 
with the numerical solution of the integral I checked that for min(gi, q2, gs) > 13 the 
difference between the truncated series and the exact result for the components of 
{X)cond is smaller than 5 x 10~^. 

In the remaining cases I employed the numerical integration described above. 



3 Theoretical Predictions 

3.1 The 2D 0(iV)- invariant vector model with iV > 2 

The 2D 0(X) -invariant vector models can formally be defined by the continuum 
action 

S = J d^x d^s{x)d^s{x), (26) 

where s is a real X-component unit vector and is the coupling constant. The 
corresponding Wilson lattice action I have already given in eq. (§). The 2D 0{N)- 
invariant vector models with N > 2 are believed to be asymptotically free |T2[ . This 
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means that the model has a finite correlation length for all finite (3 and the coupling 
g"^ vanishes with decreasing length scale. The running coupling 



2 21/ 
^ " (AT - l)e(L) 

introduced in ref. |^ is chosen such that g'^ and g^ coincide at 1-loop order of 
perturbation theory. The fiow of g^ is given by the Callan-Symanzik /5-function 

= -L-^ , (28) 

where (3{u) is given for small u by 



oo 



= -u'Y.huK (29) 

1=0 

The first three coefficients of this expansion are know from perturbation theory 0] 

6o = (iV-2)/(27r) (30) 

h = (iV-2)/(27r)2 

62 = (iV-l)(iV-2)/(27r)3. 

Of course on the lattice only finite rescalings are allowed. Therefore in ref. ^ a 
continuum step-scaling function a{s,u) is introduced, where s = L'/L is the scaling 
factor. In addition the step-scaling function on the lattice T,{s,u,a/L) is affected 
by finite lattice-spacing artefacts. The continuum step scaling function is recovered 
in the limit of vanishing a. 

\imT,(s,u,a/L) = cr(s,u) (31) 

From the analysis of the cutoff dependence of Feynman diagrams |^ one knows 
1^ that up to logarithmic corrections the leading finite lattice-spacing artefacts are 
O(a^). The dependence of the relative error 

Siu, a/L) = ^(2>^>«/-^)-^(2,^) ^ ^^^^^ ^/^^^ ^ ^^^^^ a/L)u\.. (32) 
(t(2, u) 

can be computed in perturbation theory. Some of the results given in ref. @] for 
iV = 3 are summarized in table ^ 
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3.2 The 2D XY model 



In contrast to the 2D 0(A^)-invariant vector models with N > 2 the 2D XY model 
undergoes a phase transition at a finite (3. The most accurate value for the critical 
point is Pc = 1.1197(5) []1T|. The phase transition of the 2D XY model is described 
by the Kosterlitz-Thouless (KT) theory |10|. 

For (3 < the correlation length ^i^^ik is finite and as /5 — /3c it diverges like 

6„;fcOcexp(C(/3,-/3)-i/') (33) 

where C and (3c are model dependent. The finite ^huik implies that g'^ increases 
linearly with L in the limit of large L. 

For P > (3c the theory is massless, i.e., the correlation function decays alge- 
braically in the infinite volume. The large distance behaviour in the massless phase 
is described by an effective Gaussian theory 

Seff = \(^eff I d'x{V4>{x) f , (34) 

where is the angle of the spin. The critical point is characterized by (3eff = 2/tt. 
For the Gaussian action one obtains g"^ = 1/Peff independent of L. 

The improvement of course depends on the distance r that is considered. In 
general the statistical error will increase exponentially with r. As an example I 
studied the correlation function G(r) and the effective correlation length ieffij) at 
(3 = 1.8932 and L = 64 in more detail. 

In figure 1 the relative statistical error of the correlation function is plotted as a 
function of the distance r for the three estimators. The error of the new improved 
estimator is smaller than that of the other estimators for all distances considered 
and stays almost constant in the studied region. In figure 2 the relative statistical 
error of Ce//(T) is plotted. As for the correlation function itself the error of the 
improved estimator is strictly smaller than that of the other two estimators and 
stays almost constant with increasing r, while the error of the other two estimators 
rapidly increases with increasing distance r. 

Next I tried to reproduce two of the results for the step-scaling function 

a(2, 0.7383) = 0.8166(9) (35) 

and 

a(2, 1.0595) = 1.2641(20) (36) 

given in ref. |Q. Therefore I determined S(2, 0.7383, a/L) and S(2, 1.0595, a/L) for 
various L/a in the range L/a = 6 up to L/a = 32. First one has to find for each 
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L/a the j3 at which g"^ = 0.7383 or g"^ = 1.0595 is reached. For L/a < 16 I took the 
/9's found in ref. |^ as first guess. In order to correct the small deviations from the 
wanted value of g^ 1 determined the slope of g'^{l3) from runs with smaller statistics 
(10000 measurements) at /5- values slightly smaller and larger. For the L not studied 
in ref. I started with a guess obtained from extrapolating the data for smaller 
L. Then I computed a second guess from two simulations with 10000 measurements 
in the neighbourhood of the first guess. The result of the simulation with 100000 
measurements at this second guess for f3 I give in table 1. The final correction is 
then done using the result for the slope of g'^{P). 

Next I simulated on a lattice of size L'/a = sL/a with s = 2 at the /5-values 
obtained above. The results of these simulations are also given in table 1. 

The resulting g"^ gives the finite lattice-spacing approximation J]{s,u, L/a) of 
a{s,u). In order to obtain the continuum result the limit a — > oo has to be taken. 
For 1.0595, L/a) the finite lattice-spacing artefacts even have a different sign 
than predicted by perturbation theory [Q. Nevertheless I assume, following ref. |^ 
that the finite lattice spacing artefacts vanish approximately like {a/Vf. Hence one 
might improve the convergence towards (t(s, m) with decreasing a/L by fitting the 
data for T,{s,u, L/a) according to an ansatz 

S(2, u, a/L) = (t(2, u) + Ca (a/L)^ + C4 (a/L)^... (37) 

Results for such fits are given in table 3. The statistical error of the /3's is taken 
into account in the statistical error of T,{2,u,a/L). Even when including all data 
points the /dof is small. The result for a{2, 1.0595) stays stable when data-points 
with small L are discarded; also C2 stays stable. My final estimate ct(2, 1.0595) = 
1.2639(10), taken from the fit where the two smallest L/a-values are discarded, is 
consistent with that of ref. [Q, while the statistical error is reduced by more than 
a factor of 2. Including larger lattice sizes also reduces the danger of systematical 
errors. A similar analysis leads to (t(2, 0.7383) = 0.8156(4) which is again consistent 
with ref. [g. 

Next I simulated the model at /? = 5.0 on lattices of size L = 6, 8, 12, 16, 24, 32 
and 48. The statistical error of the correlation length determined from the conven- 
tional and the cluster estimator were approximately the same. But with the new 
estimator the statistical error was reduced by a factor of about 30 compared to the 
other two estimators. The results for the running coupling are summarized in table 
4. For comparison I computed to 1, 2 and 3-loop order perturbation theory using 
the result for L=48 as start value for the integration of eq. The numerical results 
for L = 12, 16, 24 and 32 are reproduced within two standard deviations. The larger 
deviations for L = Q and 8 are well described by corrections to scaling computed in 
perturbation theory [|]. 
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Finally I also checked whether the new estimator also provides an improvement 
for the measurement of the bulk correlation length C,buik- Therefore I simulated the 
model at /3 = 1.5. The correlation length is C,buik = 11.09(2) 0]. Hence the result 
from a lattice of size L = 64 should give ^buik to a good approximation. The effective 
correlation length for the three different estimators at the distances r = 10, 20 and 
30 are given in table 5. The cluster-improved estimator is superior to the new 
improved estimator but there is still an advantage for the new improved estimator 
compared to the conventional one by a factor of 1.4 at r = 10 that grows to 2.8 at 
r = 30 in the statistical error. 



3.3 The 2D XY model. 

I simulated the 2D XY model at /3 = 0.92 in the massive phase, ai (3 = 1.1197, 



which is the best estimate for the critical coupling ITJ] and at /5 = 1.2 and 1.3 in the 
massless phase. All simulations consist of 25000 measurements. The results for the 
correlation length obtained from the different estimators are summarized in table 6. 

Let me again first discuss the reduction of the statistical error by the new im- 
proved estimator. At /? = 1.3 the statistical error of the new estimator is about 11 
times smaller than that of the conventional one and about 9 times smaller than that 
of the cluster estimator. At f3 = 1.2 the gain in the statistical accuracy is about a 
factor of 9 compared to the conventional estimator and about 6 compared to the 
cluster estimator for lattices larger than L = 8. At Pc the gain is reduced to a factor 
of 6 and 4 respectively. At j3 = 0.92 the improvement depends very much on the 
lattice size. While for the smallest lattice size L = 8 there is still an improvement 
of a factor of 2 compared to the cluster estimator, at L = 64 the cluster estimator 
is clearly better by a factor of almost 4. On the other hand the new estimator stays 
more accurate than the conventional estimator for all lattice sizes considered. 

The dependence of the improvement on the distance r is qualitatively the same 
as for the 0(3) model. Therefore I omit a detailed discussion here. 

Now let us turn to the physical interpretation of the results. In the last column 
of table 6 I give the running coupling g"^ computed from the best estimate for ^. For 
/3 = 1.2 and 1.3 the values stabilize for L > 24 within the small error bars. The 
small deviations for the smaller L are consistent with corrections that die out with a 
negative power of L. This behaviour nicely reveals the scale invariance of a massless 
theory. 

For (3 = (3c there is a tiny increase of the coupling over the whole range of 
L and the value for L = 64 is still far away from the predicted L = oo limit 7r/2. 
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However this behaviour is completely explained by 

^ - 2 " 4 1 + ^o^lnL ' ^^^^ 

which can be derived from the KT-theory Results of fits according eq. ( pHf ) 
are given in table 7. When lattices of size L < 8 are discarded the fits have an 
acceptable /dof and Qq remains stable when further data-points are discarded. 

At /3 = 0.92 g'^ clearly increases with increasing lattice size L. In order to 
measure ^}yaik one simulates lattices with L » ^buik- Here smaller distances than 
T = L are sufficient to get a good estimate for from ^e//- Therefore I evaluated 
C,eff for (3 = 0.92 and L = 64 at various distances. The results are given in table 8. 
At r = 10 the statistical error of the cluster estimator is much smaller than that of 
the new improved estimator. But still the conventional one is worse by a factor of 
1.3. 



4 Conclusion 

I showed that the new improved estimator drastically reduces the statistical error of 
the correlation length of finite systems. At values of the running coupling g'^ of the 
2D 0(3) model, such as used by Liischer et al. the error was reduced by a factor 
of 4 up to 8 compared to the cluster-improved estimator. At ^ 1/4 a reduction 
of the error by a factor of 30 is reached. 

On the other hand the error reduction seems to vanish in the limit L/S,buik 00. 
Nevertheless at L G^buik, which might be sufficient for the determination of the 
bulk correlation length, there is still an improvement of a factor of 1.3 up to 1.4 
compared to the conventional estimator, which might be interesting for models where 
no efficient cluster algorithm is known. A direct application of this method to 
gauge theories is not possible since there is no exact global symmetry that could be 
employed for the definition of an improved observable. A possible extension of the 
idea is an implementation of a multigrid scheme that has blocks with length one in 
time direction and 1,2,4 ... L in spatial direction. 

As predicted, the running coupling g'^ nicely signalled the presence of the massless 
phase of the 2D XY model. At fixed j3 > Pc the coupling g"^ rapidly approached a 
fixed point with increasing lattice size. 

In contrast to that situation, at /5 = 5.0 the g"^ of the 0(3) model decreased 
with decreasing L. The change of g^ with doubling the lattice size was more than 
100 times the statistical error. The change of g^ with varying length scale is almost 
perfectly described by the Callan-Symanzik /3-function computed to three-loop order 
in perturbation theory. 
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Figure captions 



Figure 1: The relative error of the time-shce correlation function {G{t)) com- 
puted from the conventional, the cluster-improved and the new improved estimator 
as a function of the distance r for L = 64 and P — 1.8932. 



Figure 2: The relative error of the effective correlation length computed from 
the conventional, the cluster-improved and the new improved estimator of the cor- 
relation function as a function of the distance r for L = 64 and /3 = 1.8932. 
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Table 1: Results for the correlation length ^ of the 2D 0(3) model 
obtained from various estimators of the correlation function. The sub- 
scripts con, clu and im-p denote the conventional, the cluster-improved 
and the new improved estimator respectively. In addition I give the re- 
sult for the running coupling cp' obtained from the most accurate result 
for ^. The simulations that have a direct counter-part in table 2 of ref. 
4 are marked by an x. 
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L/a 


5i 


^2 


6 


-0.00657 


0.00289 


8 


-0.00364 


0.00176 


12 


-0.00160 


0.00085 


16 


-0.00090 


0.00050 



Table 2: Perturbative results for the finite lattice spacing artefacts 
taken from ref. ||^. 



disc 


a(2, 1.0595) 


C2 


C4 


x'/dof 





1.2641(6) 


1.89(13) 


-34.7(3.0) 


0.38 


1 


1.2638(8) 


2.03(24) 


-44.(14.) 


0.34 


2 


1.2639(10) 


1.95(47) 


-37.(42.) 


0.49 



Table 3: Results of fits for S(2, 1.0595, a/ L) according to eq. 38. disc 
gives the number of discarded data points with small L. 



L 




1-loop 


2-loop 


3-loop 




6 


0.22399(4) 


0.22435 


0.22373 


0.22369 


-0.00030 


8 


0.22624(4) 


0.22667 


0.22613 


0.22609 


-0.00017 


12 


0.22959(4) 


0.23004 


0.22960 


0.22957 


-0.00007 


16 


0.23210(3) 


0.23249 


0.23213 


0.23211 


-0.00004 


24 


0.23570(3) 


0.23603 


0.23580 


0.23579 




32 


0.23836(3) 


0.23861 


0.23847 


0.23846 




48 


0.24234(3) 











Table 4: Results for the running coupling (f' aX (3 = 5.0 on lattices of 
the size L = 6 up to 48. g'^ is computed form the new improved esti- 
mator of the correlation function. The Monte Carlo result is compared 
with the result of the Callan-Symanzik /5-function to 1, 2 and 3-loop 
order. As starting point of the integration I took the MC result for 
L = 48. In the last column I give the finite lattice spacing artefact for 
a step with scaling factor 2 as computed from perturbation theory. 
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t 


Ccon 


^clu 




11 


11.023(33) 


11.031(15) 


11.031(24) 


22 


10.923(85) 


11.022(16) 


11.007(48) 


33 


10.70(26) 


11.007(21) 


11.034(94) 



Table 5: Results for the effective correlation length for various dis- 
tances t at L = 64 at /9 = 1.5. The results can be compared with 
^ = 11.09(2) obtained in ref. [| on a 128^ lattice. 



17 



fi 
p 


r 

ij 


c 

%con 


c , 


c. 


y 


l.o 


A 


7 AQf'^\ 
1 .fty ^^0 J 


7 4fi^'4'l 


7 47fi/'fi^ 


1 0701 (<X\ 
L.y) i Ui ^^y j 


L.o 





ii.4iy(^ / ) 


ii.44^0 j 


LL.oooyo ) 




L.O 


Q 

O 


iO.ZD(^iU j 


iO.io^ ( j 


iO.ZDOl^iU j 


i.U4oi^ ( j 


i.O 


1 9 
iZ 


99 Q7/'1 fi'\ 
zz.y / ^^iD ) 


99 0^(' 1 9"! 


zz. yoz^^io j 


1 ^/1f^7l'7^ 


L.O 


iO 


oU.DO^^Zi j 


"id st^n fi^ 


"in fifii n s^ 

oU.OOi (^io j 


1 n/i'?7^'fi^ 


±.»J 


9/1 
Z4 


Aft nnfrt-i \ 


/Ifi 1 


4fi nS'^('9fi'l 


1 n/iifirfi'\ 


i.O 


■^9 
OZ 


ui.ZO^oy J 


fil 18('9S'l 
Di.io^Zo j 


fil 49n(''^i^'i 


1 n49ni'fi'\ 


L.o 


4:0 


09 9'i('fi^'\ 


01 xnKn'i 


09 1 09^40 


1 041 '\((\\ 


1 


A 
ft 


fi c;^^'/!^ 


fi ^'s'i 

D.04t^^O j 


fi t^C\0(7\ 


i.zou4t^^io j 




fi 

D 


TAfP.^ 
y . 1 ^i^Dj 


7^r4'\ 
y. 1 oi^^j 


S^fil'O'l 


1 91 7^('1 1 "1 

L.Zil i Olii I 


1 9 


Q 
O 


io.iyi^y ) 


io. iy l^D j 


io.ioyi^iz J 


1 91 '^1 M 1 
i.Zioil^ii j 


1 9 
i.Z 


1 9 
iZ 


ZU.U0\^i4t j 


iy .yo^^iu J 


1 H'iOi( 1 fi^l 

iy.ooy (^io j 


1 9n07('10'l 
i.ZUy ( yLyj j 


1 9 


1 fi 
iu 


9fi zL^n 7^ 


9fi <^fi(^ 9^ 

ZO.OO^^iZ j 


9fi 41 Q^'99^ 


1911 o(^ 

i.ZiiZl iU ) 


1 9 
I.Z 


9zL 

Z4; 


zLn 1 fi('97'l 


■^O ^Ol' 1 S"! 

oy.oyi^io ) 


"iQ 74S('9S'l 
oy. / ^oi^zo J 


1 9n7fil'0'\ 

i.zu 1 Di^y j 


1 9 
i.Z 


■^9 
oZ 


c:q nO/''?1 ^ 

oo.uy\^oi j 


c;q n7l'09^ 
Oo.U ( l^ZZ J 


oz.y / o(^o4 J 


1 9081 l'S'\ 
i.ZUoil^o j 


1 9 
I.Z 




70 

i y.oo^oz j 


70 i^oMn'i 
/ y .oy J 


70 4fiSl' ^9"! 


1 9nsni's'i 

i.ZUoUl^o j 


111 07 


/I 


^ 7^(A\ 


^ 7/1 


o.oyoi^ 1 J 


i.4U4 ( (^i ( j 


111 Q7 

i.iiy ( 


D 


o.o 1 yO) 




Q KKA(^^\\ 

O.0041^iU j 


1 /in9Qi'ifi'i 
i.4;Uzyi^iD j 


111 Q7 


Q 
O 


ii.oy 1 J 


LL.ooyo J 


ii. oool^io J 


1 41 1 fi~l 

i.ftiiOliU I 


1 1 1 Q7 


1 9 

iZ 


1 fi 7fi/'i n^i 

iU. / U\^iU J 


1 fi 87l'7'l 


Ifi 01 on 0~1 
iu.y iy I iy ) 


1 41 Rf^n fi^l 

i.4;iOO\^iU } 


1 1 1 07 

L.LLu t 


1 fi 
iO 


ZZ. ftO^^iO j 


99 ^9l'0'\ 

zz.ozi^y j 


99 AAP.(9'i\ 


i.ftZOOl^iO ) 


111 07 


9/1 
Z4 


qq /I0('9'?'\ 


qq QQ/^ A\ 


qq qoni''^R^ 


1 /I '^7K(' 1 


1 1 1 Q7 


■^9 
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44 41 (^ R) 


AA '^fiqi'49'l 


1 AA')A(^A^ 


1.1197 


48 


66.33(44) 


65.92(30) 


66.082(61) 


1.4527(13) 


1.1197 


64 


88.59(59) 


87.39(35) 


87.898(89) 


1.4562(15) 


Pc 


oo 








1.570796 


0.92 


4 


3.73(3) 


3.706(12) 


3.688(6) 


2.169(4) 


0.92 


8 


6.00(6) 


6.044(18) 


6.018(11) 


2.659(3) 


0.92 


16 


8.46(9) 


8.578(25) 


8.630(25) 


3.708(11) 


0.92 


32 


10.78(28) 


10.193(26) 


10.221(86) 


6.279(16) 


0.92 


64 




10.559(42) 


8.8(1.2) 


12.12(5) 



Table 6: Results for the correlation length ^ of the 2D XY model 
obtained from different estimators of the correlation function. The sub- 
scripts con, clu and imp denote the conventional, the cluster-improved 
and the new improved estimator respectively. In addition I give the re- 
sult for the running coupling obtained from the most accurate result 
for e. 



18 



disc 





1 


2 


3 


4 




0.3466(4) 


0.3618(4) 


0.367(5) 


0.372(6) 


0.372(7) 


xVdof 


8.4 


1.4 


0.9 


0.2 


0.3 



Table 7: Results for fits of the running coupling g"^ of the 2D XY 
model at f3c following eq. 39. disc gives the number of data-points with 
small L that are discarded. 



T 


^con 






10 


10.568(24) 


10.552(8) 


10.568(18) 


20 


10.56(6) 


10.556(11) 


10.53(4) 


30 


10.66(17) 


10.537(14) 


10.55(9) 



Table 8: Results for the correlation length of the 2D XY model at 
P = 0.92 and L = 64. The effective correlation length is computed from 
the conventional con, the new improved imp and the cluster-improved 
clu estimator at the distances r = 10, 20 and 30. The results can be 
compared with ^ = 10.69(8) obtained form a 128^ lattice [p!3[ . 
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